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Abstract. We give a classification into conjugacy classes of subalgebras of 
the symmetry algebra generated by the Zabolotskaya-Khokhlov equation, and 
obtain all similarity reductions of this equation into (1 + l)-dimensional equa- 
tions. We thus show that Lie classical reduction approach may also give rise to 
more general reduced equations as those expected from the direct method of 
Clarkson and Kruskal. By transforming the determining system for the similar- 
ity variables into the equivalent adjoint system of total differential equations, 
similarity reductions to ODEs which are independent of the three arbitrary 
functions defining the symmetries are also obtained. These results are again 
compared with those obtained by the direct method of Clarkson and Kruskal, 
by finding in particular equivalence transformations mapping some of the re- 
duced equations to each other. Various families of new exact solutions are also 
derived. 



1. Introduction 

The Zabolotskaya-Khokhlov (ZK) equation is a nonlinear model of sound wave 
propagation derived from the incompressible Navier-Stokes equation [IJEllS]. The 
(2 -|- 1)- dimensional version of this equation has the form 

A{t,x,y,u) = Uxt - iuux)x - Uyy ^ 0, (1) 
and it has been studied from the Lie group approach in a number of papers 
[H HI [SI Ej • Chowdhury and Naser [1] attempted the determination of the sym- 
metry algebra of this equation and calculated some of their conservation laws. 
However, Schwarz [5] and Hereman [6] were certainly the first to provide indepen- 
dently the correct generators of the Lie symmetry algebra for this popular acoustics 
model. Although the calculation of symmetry algebras for almost all systems of 
differential equations has long been reduced to a mere function on a number of 
modern computing systems, many symmetry properties of this equation are still to 
be uncovered. 

In [7], some similarity reductions of Equation ^ to (1 -t- l)-dimensional models 
were obtained, based on the direct method of Clarkson and Kruskal 8'. More 
specifically, it was shown that if a similarity solution of Equation ^ of the form 

u = U{t,x,y,W{^,r])), £,=^{t,x,y), 'n = 7j{t,x,y) (2) 

can be found by solving a reduced (1 + l)-dimensional equation, then, when 7^ 0, 
such an equation must be one of three non-equivalent (1 -I- l)-dimensional equations 
found in that paper. However, these reduced equations are either determined only 
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implicitly in terms of solutions of certain partial differential equations (pdes), or 
they depend on up to four arbitrary functions, and the same can be said about 
their solutions. In the same paper, using again the direct method of Clarkson and 
Kruskal and a restrictive anzatz, the most general ode that every similarity solution 
of II]) obtained by solving an ode must satisfy is shown to be of the form 



w'"^ + ww" + {Az + B)w + 2A^ -{Az + Bf. (3) 

Some of the unanswered questions raised in [7] were how to find all the non equiv- 
alent similarity reductions of the ZK equation to an ode by the classical Lie sym- 
metry approach, and whether there is any connection between these two types of 
reduction techniques. The first of these two questions stems from the fact that 
the symmetries of Equation ([1]) depend on three arbitrary functions of time, and 
so its similarity reductions are usually achieved by restricting these functions to 
some specific types of elementary functions, such as exponential or simple polyno- 
mial functions of time . It also stems from the fact that no classification of 
low-dimcnsionial subalgebras of the ZK symmetry algebra into conjugacy classes is 
available. 

In this paper, we obtain all canonical forms of non-equivalent one- and two- 
dimensional subalgebras of the symmetry algebra L of ([T]), under the adjoint rep- 
resentation of the symmetry group. We then apply them to obtain all similarity 
reductions of the ZK equation to (1 -I- l)-dimensional equations, using the classi- 
cal Lie approach. The same reductions are also obtained by direct case analysis. 
We thus show that in addition to its simpler algorithm and other properties. Lie 
classical reduction method gives rise not only to simpler equations, but it may also 
yield more unified and general reduced equations than those expected from the di- 
rect method (of Clarkson Kruskal). Next, by transforming the determining system 
for the similarity variables into the equivalent adjoint system of total differential 
equations, similarity reductions to large classes of odes which are independent of 
the three arbitrary functions defining the symmetries are also obtained. The lat- 
ter system of total differential equations allows for an easier determination of the 
invariant functions defining the similarity coordinates. Large classes of similarity 
solutions depending on much less arbitrary functions than those obtained in [7] are 
also derived in this way. 

Finally, we find equivalence transformations mapping some of the equations that 
we've obtained by the Lie classical method to some sub-equations of the reduced 
equation ^ obtained by the direct method of Clarkson and Kruskal. Our discus- 
sions also show that in principle any reduced equation achievable with the direct 
method can also be achieved by Lie classical method although the converse is totally 
out of question, as far as the properties of the reduced equations are concerned. 

This paper is organized as follows. In the next section we discuss the symmetry 
algebra of the ZK equation and determine its algebraic structure as well as its 
connection with Kac-Moody-Virasoro (kmv) algebras. Section [3] is devoted to the 
classification of low- dimensional subalgebras of L and Section |4] to the similarity 
reductions of the ZK equation. We investigate the connection between the two 
types of reductions invoked above in Section [5] . Some concluding remarks are 
given in the last section. 
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2. Symmetry group of the ZK equation 

2.1. Structure of the symmetry algebra. The Lie algebra of the ZK equation 
is well known [SI [S] . This is the Lie algebra defined by the infinitesimal generators 
of the point symmetry group G of the equation, that is, the Lie group of point 
transformations that map every solution of the equation to another solution of the 
same equation |9l [10] . These infinitesimal generators are vector fields of the form 

V = (,i{t,x,y,u)dt + ^2{t,x,y,u)dx + ^3{t, x,y,u)dy + (p{t,x,y,u)du (4) 

acting on the space of independent variables coordinatized by {t, x, y) and the space 
of the dependent variable coordinatized by u, and such that the second prolongation 
Pr*-^-* V of V satisfies 

Pr^^K A{t,x,y,u) =0. (5) 

A{t,x,y,u)=0 

Equation ^ completely determines the unknown functions ■^1,^2,^3,0 defining 
V |9l [22l [To J . For the ZK equation, the linearly independent vector fields, i.e. the 
generators of the Lie symmetry algebra are given by the operators 

vo = 2xdx + ydy + 2udu (6a) 

Xg = gdx - g'du (6b) 



Yh = 7;yh'dx + hdy - -yh"du (6c) 



2/ ^ f9t + \ {2xf + y^f") + ^fdy + i {^4ur - 2xf" - y'f") du (6d) 

where /, g, h are arbitrary C°° functions of the time variable t defined on some open 
subset of K, assumed to be the based field of L, and where a prime represents a deriv- 
ative with respect to t. The ZK symmetry algebra is therefore infinite-dimensional, 
and its commutation relations are easily found to be as follows 

[vo,Xg] = -2xg, [xg,y/i] = (7a) 

[vo,y/i] = -Yh, [xg,z/] = X(/,g/3_/g/) (7b) 

[vo,z/] = 0, [yh,Zf] ^ y(2^.,^_y;,,) (7c) 

[Xgi,XgJ = [yfti,y/i2] = ^ih^h'^^h[h2)/2 (7d) 

[z/i,z/2] = z^f^f^_f^f^). (7e) 

These commutation relations show that the ZK symmetry algebra L has the struc- 
ture of a semi-direct sum Lie algebra, L = TZ + S , where TZ = (vo,Xg,y/i) is the 
radical and S ^ (zf) is the semisimple part, also called Levi factor of L. It is indeed 
a well-known fact |11] that the infinite dimensional Lie algebra generated by is a 
simple Lie algebra, i.e. it has no nontrivial ideal. This is easily established by first 
nothing that by a result of Cartan [H], the Lie algebra J(/) = {fit)dt : / G C°°(/)} 
of vector fields on the open interval / of R is a simple algebra, and that the map 

ct: J(/) -> 5: f{t)dt ^ Zf(t) 

is a Lie algebra isomorphism. We also note that the radical TZ, which is solvable by 
definition, is actually nonnilpotent. Its nilradical, i.e. its maximal nilpotent ideal 
is generated by (xg,y/,) . 

The commutation relations ([7]) also show that the infinite-dimensional subalgebra 
Lq — {xg,yh,Zf) of L can be identified with a subalgebra of a KMV type algebra. 
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Indeed, by restricting the arbitrary functions /, and h to be Laurent polynomials, 
the resulting commutation relations take the form 

[xt„.,xtn] = 0, [xt..,yt.] = (8a) 

[Zf,„,zt,.] = (n - m)zt,„+„-i, [xt,^,Zt,.] = x^^+^-i (8b) 

2rn — 3n n — m 

[yf™,zt„] = yt™+„_i, [y^,"^,yt-] = — - — x^^+^-i, (8c) 

and this shows that (xg,y/i) generates the corresponding Kac-Moody subalgebra 
of Lq, while (zf) generates the Virasoro subalgebra These algebraic prop- 

erties of the ZK equation, which can be linearized by a generalized hodograph 
transformation [18j . are in agreement with a widespread belief about integrable 
(2 + l)-dimensional equations. Indeed, the symmetry algebra of most integrable 
(2 + l)-dimensional equations are infinite-dimensional and can be identified with 
a subalgebra of a kmv type algebra [ill US HI] . There are nonetheless some ex- 
ceptions provided for instance by the breaking soliton equation and the Zakharov- 
Strachan equation [15] which are integrable but do not admit a kmv type symmetry 
algebra. On the other hand, the symmetry algebra of every known non-integrable 
equations has no kmv structure. This suggests that kmv structures are in some way 
which is still to be clarified associated with integrability. By an integrable equation 
here, we refer to equations allowing a Lax pair, an infinity of conservation laws, 
soliton and multisoliton solutions, a family of analytic periodic and quasi-periodic 
solutions, as well as a number of similar properties. 

2.2. Group transformations. One of the most important properties of the sym- 
metry group of a differential equation is to transform one solution of the equation 
into another solution of the same differential equation. These group transforma- 
tions are generated by some one-parameter group of transformations, each of which 
is the local flow -0(e,w), where w — {t,x,y,u) and e G M, of a basis element v 
of the Lie symmetry algebra, acting on the space of independent and dependent 
variables. More specifically, the local flow ?/'(e, w) of v is given for every e for which 
it is defined by 

5eV'(e,w) = ■0(e,w) = v|^^^,^^, and V(0,w)=vir, (9) 

where = d/de. The full group action can be obtained by combining all the 
one-parameter group actions determined by each infinitesimal generator. We will 
use the notation 

V'(e,w) = w = {i,x,y,u), 

and we let the generic generator of the ZK symmetry algebra have the form 

V = fcovo +Xg +y,j +z/, (10) 

where fep is a scalar. We also denote by F' the derivative of a function F = F{t) of 
the time variable t. To find the flow of v for the ZK equation, we have to distinguish 
a number of cases. 



Case 1 : g — h — f — 0, and fco 7^ 
In this case v = vq, and by assuming without loss of generality that fco = 1, we 
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have 

dei^O dey = y 

dfX = 2x deU = u. 

Consequently, the corresponding group action is given by 

•^(e, w) = {t,xe'^^,ye^,ue^'^) = w. 

It then foUows that 

u{t,x,y) = e^^u{t,xe~'^'^,ye~'^) 
is a solution of the ZK equation whenever u{t, x, y) is a solution. 

Case 2 : f = h = 0, g 0, ko is arbitrary 
We have in this case v = Aiovo + Xg. 

Case 2a : fco = 
We find that 

w = {t, x + eg, y,u- eg'). 

Consequently, 

u{t, X, y) = u{t, X - eg, y) - eg' 
is a solution whenever u{t, x, y) is. 

Case 2b : fco 7^ 
In this case, w is given by 

i = t 

x=T^ie"'°'{2kox + g)-g) 
2fco 

y = ye''°' 

u=^y'^%2kou-g')+g'). 
Therefore, if u{t, x, y) is a solution, then 



u{t,x,y) = ^^[e'''''\2kou*-g')+g'] 



is also a solution, where 



u =u 



t,^^{{2kox + g)e-^^°'-g),ye-''-' 



Case 3 : f = 0,h and g and fco are arbitrary 
We have two subcases in this case. 

Case 3a : fco = 
We have 

w = {t,x + ge + e{y + he)h', y + 2he, u - e{g' + {y + he)h")). 
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Consequently, whenever u(t,x,y) is a solution, so is 

u{t, X, y) = u* - e{g' + y{-ht)h") 

where 

u* = X — ge + e(— y + he)h\ y — 2he). 
Case 3b : ko ^ 

We find that w is given by 
i = t 

i = [-kag + 2hh" + e^°\-Ahh' - 2koyh') + e^^'^'^ikog + 2klx + 2hh' + 2koyh')\ 



y=^[-2h + e^^'%2h + koy)\ 

w — ^ [-fco.g' + 2hh" + e'""'{-Ahh" - 2koyh") + e^'"''{-2k^u + k^g' + 2hh" + 2kQyh")\ 

2fcQ 

This shows that whenever u{t, x, y) is a solution, so is 
^(^~' y) = i [2fc2e2^-«^w(i, x, y) - fco(-l + e''"^).9' + 2(e^-«^ - + e"«^)/i - fcoy)/^"] 

where 

t = i 

[e-'"''^(fco(ff - e'^ng + 2fco5;) + 2/i'(-l + e'="^)((e'=«^ - - k^y)] 

y = -^ [e'''°'i-2h + 2e''°'h - koy)] . 
fco 

Finally, we are left with the case where / 7^ 0. This case leads to relatively 
complicated or too long formulas, and we omit it here. 

3. Classification of low-dimensional Symmetry algebras 

Group-invariant solutions corresponding to two subgroups which are conjugate 
under the adjoint action of the symmetry group lie on the same orbit. There- 
fore, in order to know all the similarity solutions of the ZK equation invariant 
under s-parameter symmetry subgroups, it is sufficient to have a classification of 
s-parameter subalgebras, and an expression for the adjoint representation Ad of G. 
In the case of the ZK equation, as explained in the next section, we shall only need 
a classification of s-parameter subalgebras under the corresponding adjoint action 
ad of the Lie algebra L of G on itself, where 1 < s < 3, and where ad denotes the 
differential of the Ad . By the connectedness of G, in order to find an expression 
for Ad, we only need to find it for the flow Ad(exp(£l^)) • v through v of the one- 
parameter subgroup eyip{eVi) generated by each basis element Vi of L. However, if 
we denote by 14d the infinitesimal generator of Ad corresponding to each V (z L, 
then Vad coincides with the Lie bracket on L: Vad • v — [v, V] for all v e L. We 
can therefore reconstruct the adjoint representation of the group from that of its 
Lie algebra, and formulas for Ad(exp(eF)) • v based on Lie series or on properties 
of the flow exist Page 205]. 
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For the purpose of finding a representative list of all similarity solutions of the 
ZK equation, we shall only need as indicated to classify one- and two-dimensional 
subalgebras. Techniques for performing these classifications can be found in [161 
1171 [9], They generally consist in mapping or 'reducing' a system of Lie subalgebra 
generators to an equivalent system until a canonical representative is achieved. For 
the ZK symmetry algebra, we shall need the following set of lemmas which are 
essentially based of the commutation relations 

We shall often denote collectively by Wg, where q is some function of time, all 
of the three families of generators :x-g,yh and Zf in It is clear that is a 
linear function of its argument q. When g is a constant function represented by its 
numeric value q, Wg will be denoted by W(g). We shall also often represent a general 
vector of the form (llOp simply by its components in terms of fcoi Xg, y/,, and Zf. 

V = fcovo +Xg +y,i + z/ = {fco,5,^, /} (H) 

Lemma 1. Let f,g and h be given functions oft, and ko a given constant. Denote 
by a and f3 some numbers and by G and H some functions of t. 

(a) For some a and G, vq + Xg is mapped to Vq wnrfer Ad(exp(a xq)). 

(b) For some a, /3, G and H , vq + Xg + is mapped to vq under 
Ad(exp(Q;XG)) Ad(exp(/3 yn))- 

(c) For some a, /3, G and H , /cqVo + Xg + y/i + zj is mapped to fegvo + z/ under 
Ad(exp(Q!XG)) Ad(exp(/3yif )), provided that + f"^ ^ 0. 

(d) For some (3 and H , Xg-|-y/i is mapped to either Xg ory^ under Ad{exp((3 y h)) 

(e) For some a and G, voH-Xg+y^ is mapped to VQ-t-y^ under Ad{exTp{a xg)) . 

(f ) For some a and G, ko'VQ + Xg + y/j + zj is mapped to fcoVo + y/i + z/ under 
Ad(exp(axG)), provided that k^ + f'^ 0. 

Proof. We have Ad(exp(Q;XG))(vo) = vq — 2q;xg, and Ad(exp(axG))(vo -t- Xg) = 
Vq — 2axG + Xg, which has component {1, g — 2aG,0,0}. It suffices then to 
choose G — g/{2a) with a ^ 0, which proves part (a). For part (b), we no- 
ticed that when h is zero, the statement reduces to that of part (a). Otherwise, 
V = Vq + Xg + yh is transformed under Ad(exp(/3y//)) to a vector with component 
{1, g + {l3/2){-Hh' + hH'),h - (3H, 0} . If we therefore choose H = h/(3 to cancel 
the term yh in V, the result will follow from that of part (a). For (c), we note that 
if / = 0, the problem reduces to that of part (6). Thus we first assume that / ^ 0, 
and /cq = 0. Then Ad(exp(/3y/f )) will transform = Xg + y/i + z/ into a vector 
V'^^^ of the form {O, gi, ft- -f + /i?'), /} , where gi is a certain function 
of t. Thus to cancel the term yh in V, it suffices to choose H as a. solution to the 
equation h + l3{—^H f + fH') = 0. Notice that this is a hnear first order differential 
equation of the general form 

K'{t) = Q{t,K{t)), (12) 

where Q is a linear function of the unknown function K, and such linear equations 
always have a solution. For instance, the latter equation for H has solution 

H{t) = Gfitf' + fitf' -^rf^, (/3 ^ 0), 

where C is an arbitrary constant. Now suppose that V^-^'^ reduces to a vector of the 
form {0,g1,0j} when H is so chosen that the third component of vanishes. 
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Then Ad(exp(axG)) will map the resulting vector V'^^'> to 

={0,g- a(-(l/3)G/' + /G'), 0, /} . 

The condition that the second component of F*^^^ vanishes is a linear equation in G 
of the form (fT2|) . showing that we can reduce V to x.f = fcoVg + x/. If on the other 
hand both ko and / are nonzero, then Ad(exp(/3yij)) will map V — fcovo + x^ + 
Yh + Zf to a vector T/^^^ with component {/cq, 5i, ^ + /? (--f (3 + 2/') + fH') , /} , 
for some function t/i depending on H. Let _ffo be a function for which the third 
component of V^-^"^ vanishes, and denote again by = { ^0:51,0,/} the corre- 
sponding image of V under Ad(exp(/3yHo))- Then, 

Ad(exp(axG))(l^(^)) = {ko,g - (l/3)a(G(6 + /') - 3/G'), /} , 

showing that we can map V as indicated to fcoVo + 2/, and this proves (c). 

For (d), if h — 0, then V = Xg +y/i reduces to Xg and we choose (3^0, otherwise 
it suffices to choose H as a solution to the linear equation 

g+lf3{-Hh' + hH')^0, 

in order to map V to y/j. For (e), it suffices to choose G — g/{2a), with a 7^ 0. To 
prove (f), we simply note that if kg = 0, we let G be the any solution to the linear 

ODE 

5 - a Qg/' + /G'^ = 0, {a^O). 

If / = 0, we let G be a solution to g — 2aG — 0, and if both feg and / are nonzero, 
we let G be the solution to the linear ode 

g-aQG(6 + /')-3/G'^ =0, (a^O). 
This completes the proof of Lemma [T] □ 



Lemma 2. Let g, h and f be given functions of time. Then for any e G M, and for 

every nonzero function K — K(t) , we have 



(a) 



Ad(exp(ezK))(xg) = :xiGit,e)i Ad(exp(£Zi<-))(y,,) = yH(t.,e), 
Ad(exp(ez/f))(z/) = ZF(t,e), 



(13) 



SYMMETRY PROPERTIES OF AN ACOUSTIC**(S) MODEL 



9 



where the functions G(t,e),H(t,e) and F{t,e) are given by 

G{t,e) ^ KUy/^Qi (e- [ , and Qi satisfies 

V Ji K{s)J 



(14a) 



git) = Kity^'Q, [- f (14b) 



and Q2 satisfies (14c) 



K{s) 
* ds 



K{s) 



* ds 



hit) = Kit)'/'Q2 ( - ) . (14d) 



* ds 



1 



Fit,e) = Kit)Qz \ I ^ ) ' '^"^ satisfies (14e) 



/(0 = i^(0Q3(-/^). (14f) 

(b) Lei Wq, where q = 9(^)1 denote collectively all generators of the form 
Xg,y;j and Zf. Then, whenever q ^ 0, there exists a function K such that 
Ad(exp(£Zif))w, = W(i). 

Proof. From Equation ([7]) we have [zi<:,Xg] = ^{Kg' -K' g/'i)^ and this shows that if 
we express Ad(exp(e z/f ))(xg) in terms of the Lie series and use the hnearity of Xg 
as a function of its argument g, then Ad(exp(ezx))(xg) must be of the form 'x.G{t,e) 
for a certain function G = Git,e). Now, the expression for Git,e) foUows from the 
properties of the flow Ad(exp(ezi<-))xg of Ad through Xg under the one-parameter 
subgroup exp(ez;f) generated by K. The expressions for Ad(exp(£ Zi<-))(y?i) and 
Ad(exp(ezif))(z/) and for the corresponding functions Hit,e) and Fit,e) are de- 
rived in a similar way, and this proves (a). 

Part (b) simply says that if for instance we have ^ 0, then we can choose the 
function Kit) in such a way that the resulting function G(t, e) in is 1, and that 
the same holds for both h and / and for the corresponding functions i?(t, e) and 
Fit, e) in (|13p. The existence of such functions K is just a consequence of a result 
of Neuman [19] . A proof of the existence of K is also given in [11] . This completes 
the proof of the lemma. □ 



3.1. Classification of one-dimensional subalgebras. To implement this clas- 
sification, we denote as usual by v = fegvo + Xg + y/i + z/ a general nonzero vector 
in L, and for Vi,V2 G L we write vi ^ V2 if V2 = Ads(vi), for some group el- 
ement B in G. When both and / equal zero, we have v ^ Xg or v ^ y/j by 
part (e) of Lemma [TJ Otherwise, we have + f"^ ^ 0, and hence v ^ fcgVo -|- zj 
by part (c) of Lemma [1] Thus every one-dimensional subalgebra of L is equiva- 
lent under the adjoint representation to either Xg, y^ or fcoVg + z/. Since we have 
AdB(fcoVo -I- z/) = fcoVo -I- Adsi'^'f), fo^' every group element B, it follows from 
Lemma|2]that if / = 0, then (fcoVo-l-z/) ^ vq, other wise (fcoVo-|-z/) ^ fcoVo-l-Z(i). 
The same lemma thus implies that every one dimensional subalgebra of the ZK sym- 
metry algebra is equivalent under the adjoint representation to one of the following 
Lie algebras 
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-Ci,! = {x(i)} , £i,2 = {y(i)}, 'Ci,3 = {vo} or £1,4 = {fcovo + z(i)} . 

(15) 

Note that according to ([6]), we have 

X(i) = da;, y(i) = dy, and Z(i) = dt, (16) 

which shows that the canonical forms thus obtained are very simplified. To see 
why these four Lie algebras are non equivalent, we first note that Ad(exp(£Vo)) 
acts only diagonally, by scalling its argument, while for every function F — F{t), 
Ad(exp(e z^? ) ) leaves vq invariant and maps, by Lemma[2l Xg , and zj to Xg , yp, , Zj , 
respectively, for some functions g — Q{t,e), f) = i){t,e) and f = f(i,e). On the other 
hand, if we set 

= Ad(exp(Q;xG)) Ad(exp(/3yi/)), 

then we have, 

T{vo) = {l,-2aG,-PH,0} ^ ^ro 

and 

^(fcovo + Z(i)) = {fcovo, a{-2G + G') + {l/A)f3^{H'^ - HH"),P{~H + H'), 1} 
kovo + Z(i). 

We also have J^(x(i)) — X(i), and .F(y(i)) = y(i)- This shows that the four Lie 
algebras thus obtained are non-equivalent under the adjoint representation of G, 
and thus completes the classification of one-dimensional subalgebras of L. 

3.2. Classification of two-dimensional subalgebras. Every two dimensional 
subalgebra L2 of L is solvable and if we denote by S = {Vi, V2} a basis of L2, then 
either L2 is abelian and thus [Vi, V2] = 0, or L2 has a nonzero nilradical, in which 
case its commutation relations can be put in the form [l^i, V2] = Vi. To classify two 
dimensional subalgebras of L, we let Vi be in one of the canonical forms (fT5| , while 
V2 = V is a general vector of the form pT|) . By the possible forms for Vi given 
in p^ . we shall therefore have to consider four basic cases for B. To reduce L2 to 
a canonical form, we transform it under maps of the form Ad^, for appropriately 
chosen group elements B in the symmetry group G, and also make use of the 
commutation relations in L2 to obtain some restrictions on the functions g, h, f 
and the free parameter ko defining V2 = v. If Bi is a basis of a subalgebra and 
B2 a basis of a subalgebra S2, we write Bi ^ B2 if S2 — AdB(S'i) for some B £ G. 
In this subsection, we shall let ci,C2 and C3 denote three arbitrary constants. 

3.2.1. Abelian subalgebras. 
Case (a): B = {vo,v} . 

The condition [vg, v] = 0, implies that g = h = 0. Thus v = fcgVo + zy. Since for 
every e G M and every function K we have 

Ad(exp(ezK))(fcoVo + z/) = /cqVo + Ad(exp(ezif))(z/), 

and since Ad(exp(ez/f ))vo = vq, it follows that if / = 0, then v ^ vq. Otherwise, 
V ~ A;oVo -I- Z(i-) . Since L2 is two-dimensional, we must have B ~ {vq, ^oVq -I- Z(i) } . 

Case (b): S {/ciVq -|- Z(i), v} , /ci e M. 
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We must have in this case 

g^cie^^'\ h = C2e^'\ and/==C3. (17) 
Thus B ^ {fcivo + Z(i), v} , where the components g, ft, / of v are given by (|17p . 
Case (c): S={x(i),v}. 

We must have / — Gkot + C3. If fco = and C3 = 0, then by part (d) of 
Lemma [U under Ad(exp(/3y^f )), we have v ^ Xg or v ^ y^, and thus, since 
Ad(exp(/3y//))xg = Xg for aU g, we have B ^ {x(i),Xg} or B ^ {x(i),y,i} . If 
fco 7^ or C3 7^ 0, then by part (c) of Lemma[l] we have v ~ fcovo + Z(6fcot+c3)' and 
hence 

B {x(i),A:oVo +Z(efc„4+C3)} ■ 
Case (d): S={y(i),v}. 

We must have / = — (3/2)fcot + C3, and ft = C2. Thus if fco = and C3 ~ 0, 
then V ^ Xg +y(c2)j and no further reduction of the basis B = {y(i),Xg + y(c2)} is 
possible. If fco 7^ or C3 ^ 0, then 

S ~ {y(l), fcoVo + y(c2) + Z(-(3/2)fcot+C3)} • 

This is because by part (f ) of Lemma [l] we have v ~ fcovo + yc2 + ^/ j under 
Ad(exp(a Xg)), and the operator Ad(exp(axG)) leaves y(i) invariant for every a 
and G. Again, all the canonical representatives of two-dimensional subalgebras thus 
obtained are clearly pairwise nonequivalent, and this completes the classification 
problem in the abelian case. 

3.2.2. Non-abelian subalgebras. 

Case (a): S = {vo, v} or S = {fcivo + Z(i), v} , /ci € R. 

Since Vq is not in the derived subalgebra of L, the case B ~ {vo,v} cannot 
occur, while the case B — {fciVp + Z(i), v} occurs only if ki — 0, in which case we 
must have 

g — ci, h — C2, and f = t + 03 
Consequently, in this case 

B - {z(i), fcovo + X(c^) + y(c^) + Z(f+C3)} . 
Case (b): S={x(i),v}. 

We must have in this case / = (3 — 6fco)i + C3. Thus if fco = 1/2 and C3 — 0, 
then by Lemma[Tl part (b), Ad(exp(axG)) Ad(exp(/3y//)) will map v to vq and 
leave X(i) unchanged. Hence, B ^ {x(i),vo} in this case. If either kg ^ 1/2 
or C3 ^ 0, the same lemma implies that v ^ fcoVo + zj, under a similar trans- 
formation that leaves X(i) unchanged. Consequently, in this second case we have 

13 ^ {X(i), /coVo + Z((3„6feo)t+C3)} • 

Case (c): S={y(i),v}. 

The commutation relations imply that / = (3/2)(l — ko)t + C3, and ft = C2. Thus if 
fco = 1 and C3 = 0, by part (e) of Lemma [U Ad(exp(axG)) maps v to vo + y(c2)- 
Since this operator leaves y(i) invariant, we must have B ~ {y(i),vo +y(c2)} in 
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this case. If fco 7^ 1 or C3 7^ 0, part (f) of the same lemma shows that Ad(exp(a xq)) 
maps V to fcovo + y(c2) + ^Z- Consequently, 

B {y{l), fcoVo + y(c2) + Z((3/2)(i_fco)t+C3) } ■ 

Since the canonical forms thus obtained are non-equivalent by construction, this 
completes the classification of two-dimensional subalgebras of L in the non-abelian 
case. 

We have thus obtained the following list of canonical forms of non-equivalent 
two-dimensional subalgebras of L, in which fco, fci, ci, C2, and C3 are free parameters, 
unless otherwise specified. 



(1) Abelian subalgebras 



-^2,1 


{vo, fcoVo + Z(i)} 






^2,2 


= {/civo -f Z(i),A:oVo 4 




-=2(C3)} 


-^2,3 


= {x(i),x<,}, (5V 0) 




-^2,4 


- {x(i),y^} 






-^2,5 


= {x(i), fcoVo -1- Z(6feotH 






•^2,6 








^2,7 


= {y(i), fcovo + y(c2) H 


- Z((_3/2)(feot+c3))} , 




(2) Non 


— abelian subalgebras 






•^2,8 


= {z(i), fcoVo + X(ci) 4 


-y(c2) +Z(t + C3)} 




^2,9 


= {x(i),vo} 






-^2,10 


= {X(l) , fcoVo + Z((3_gfeo)(+c3) } , (^0-5) 




-^2,11 


= {y(i),vo +y(c2)} 






^2,12 


= {y(i),fcovo + y(c2) H 


- Z((3/2)(i-fco)t+c3)} , 


iko-ir+4 



4. Similarity reductions of the ZK equation 

We turn our attention in this section to the problem of finding the group- invariant 
solutions of the ZK equation. By group invariant solutions we mean solutions which 
are invariant by the group transformations in the sense that, roughly speaking, each 
of them is transformed into itself by every group transformation. In other words, 
their graph is a locally G-invariant subset. In the case of a PDE such solutions 
can be found by solving a differential equation in fewer independent variables. We 
shall therefore seek a reduction of the ZK equation to (1 -I- l)-dimensional equations 
or to ordinary differential equations. To reduce the equation to one with s fewer 
independent variables, the general procedure [3] is to look for subgroups of the full 
symmetry group whose orbits have dimension s. Each such subgroup yields a set of 
3 — s functionally independent invariants which can be written in the form ^ = r7(w), 
where w — (t,x,y), in addition to another functionally independent invariant of 
the form w — C(w,m). We must clearly have 1 < s < 3, since the ZK equation has 
exactly three independent variables. Solving this last equality for the solution u of 
the original equation shows that u will always be determined in the case of the ZK 
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equation by some equations of the form 

u{t,x,y)^a + /3w{0, C = f?(w) (18) 

where a and /3 are functions of w. Equation IjlSp can be used to rewrite the original 
equation in terms of the invariant functions ^ = r]{w), and w considered as new 
variables, and this yields the reduced system in only 3 — s independent variables 
Expressing the solution of the reduced equation in the form w = F{^), and 
substituting in this last equality the expressions for w and ^ from equation 
the solution u of the original equation can be found. We treat separately reductions 
by one dimensional subgroups and reductions by two dimensional subgroups. 

4.1. Reduction by one-dimensional subgroups. As already indicated, a re- 
duction by one dimensional subgroups will yield a partial differential equation in 
only 2 = 3—1 independent variables. The generic form of the infinitesimal genera- 
tor of such a group is given by (|10|) . In the case of a reduction by one-dimensional 
subgroups, the function the ^ = ?7(w) in (|18p will have the form ^ = (771,772), and 
for simplicity of notation we shall set r ^ rji and z = 772. In order to completely 
define a reduced equation, wc will only need to give the explicit formulas for u, z, 
and r corresponding to Equation and the reduced equation itself. We first give 
reductions based on the classification of onc-dimcnsional subalgebras given in (fTS]) . 
and then reductions based of a direct case analysis afterward. 

4.1.1. Reductions based on one- dimensional subalgebras classification. This proce- 
dure consists in finding the reduced equation corresponding to each of the four 
non-equivalent canonical forms of one-dimensional subalgebras of L given in ()15|1 . 
The list of all solutions to the four resulting reduced equations should represent an 
optimal list of 1-parameter group invariant solutions, with the property that every 
other 1-parameter group invariant solution of the ZK equation can be mapped to 
precisely one solution in the list via the adjoint representation of G. 

(1) Reduction by £11 = {x(i)} . 

For X(i) = dx, the reduction formula is simply u — u{t, y), r ~ t, and z = y. 
The reduced equation is the linear equation 

with solution u = qiy + q2, where qi and (72 are arbitrary functions of t. 

(2) Reduction by £i_2 = {y(i)} • 

For y(i') = dy, the reduction formula is u = u{t,x), r = t, z = x. The 
reduced equation is 

~ul + Ut^x - UUx,x = 0. 

(3) Reduction by £1^3 = {vq} . 
The reduction formula is 

u — xw{t,z), r = t, z = y^/2, (19) 
and the corresponding reduced equation is 

Wr — 2Wz — {w — ZWz)^ — ZWr,z — ^ZWz.z — Z^WW^.z — 0. (20) 
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(21) 



(4) Reduction by £1,4 = {/cqvo + . 
The reduction formula is 

u = xw{r,z), r — 2kot — \n{x) , z = y'^/x, 
which gives rise to the reduced equation 

2Wz + { — W + Wr + ZWzf' + 2fco(Wr,r + ZWr,z) + ^ZWz,z 
+ w{ — Wr + Wr,r + z{2Wr,z + ZWz.z)) — 2k()Wr = 0. 

For fco — 0; this last equation corresponds to the much simpler reduction 
by {z(i)} given by 

ul + uux^x + Uy^y = 0. (22) 

4.1.2. Reduction by direct case analysis. This procedure consists in obtaining all 
possible reductions of the ZK equation by one-parameter subgroups, by a suitable 
consideration of each relevant case separately. Although it gives directly all possible 
cases of reduced equations, it might involve tedious or too long calculations in 
certain cases. 

Case 1 : g = h = f = Q, and fco 7^ 
This is just a reduction by the subgroup generated by vq. The reduction formula 
and corresponding reduced equations are already given in (jl9p and ()20|) . 

Case 2 : h ~ f ~ 0, and 5 7^ 0, and fco arbitrary 

Case 2a : fco = 
In this case we have 

u——{w — x), r = t, z = y. (23) 
9 

The reduced equation is 

g" + g'wz,z = 0. (24) 

This gives rise to the invariant solution 

u = —{w — x), w = a{t)y + 6(t), g{t) — ait + a2 
9 

where a{t), b{t) are arbitrary functions, while ai and 02 are arbitrary constants. 

Case 2b : fco 7^ 
The reduction formula is 

u^y^w{t,z) + ^, r = t, z = \{x + g/i2ko)), (25) 

and the reduced equation is 

~2w + 2zwz ~wl + Wz,t - (w + 4:z'^)wz,z = 0. (26) 
Case 3 : / = 0, and 7^ 

Case 3a : fco = 
The reduction formula is 

u=^{w{t,z)~{2yg' + y'h")), r = t, z = x - ^{g + yh' /2) . (27) 
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(28) 



This system implies that 



w 



h — —t + aia2t + as, 



-t + ai 

and the corresponding invariant solution is 

-t{2ai + t)xa2 - 2xa3 + a{t){ai + t){t{2ai + t)a2 + 2a3) + y{g + (-oi - t)g') 

u - 

(ai +t)(<(2ai +i)a2 + 2a3) 

where a(t) is an arbitrary function, while oi, 02 and 03 are arbitrary constants. 

Case 3b : fco 7^ 
The reduction formula is 



u^[2h + kQyfw{t,z) + 
1 



k^g' + 2{h + kay)h" 



2k1 



kog + 2{h + koy)h' 



2kQ 



{2h + koyf 

and the corresponding reduced equation in which we may set k^ = 1 in 

- 2klw + 2klzwz -wl+ Wz,t - {w + Aklz'^)wz,z = 0. (29) 
Case 4 : / 7^ 

In this case the coefficient of dt in the expression of the operator z/ is nonzero 
and thus we may no longer consider the variable t as trivially invariant. We make 
the simplifying assumption that in (jlOp we have k^ = 0, so that we are looking for 
a reduction of the equation by a generic vector field of the form 



Yh 



(30) 



The expressions for the invariant functions w.z, and r are found after some long 
calculations to be given by 



1 



54/4/3 
1 

+ 54 



-9 



54x/ + 3/2/3/iG - 9y{h + yf 
g Gh' 1 



G^f" \ dt + G'^f 



2 rl 



fi/3 3J2/3 9 

367./^ -(h + 2yf'f + 6/(.g + 2xf' + y{h' + yf")) 
36/4/3 



(31a) 



(31b) 



(31c) 
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where 

The reduced equation in this case is 

wl + WWz,z + Wr. r = 0. (32) 

It appears that equations and ([5^ are exactly the same, despite the long 
calculations leading to the latter equation, and this confirms in some sense the 
results of part (c) of Lemma[l]and that of Lemma[5]which assert that when fco = 0, 
and / 7^ 0, the corresponding vector v of (|30p is equivalent to Z(i)- In a similar way, 
in the case where / 7^ and fco 7^ 0, the equation obtained by direct case analysis, 
although not calculated in this paper, should match that obtained in (|2ip . Similar 
correspondences can be established between the equations obtained in Section [4. 1.1 1 
and Section l4X2l 

The classification obtained by direct case analysis in this section contains an 
exhaustive classification of all possible similarity reductions by one dimensional 
subgroups, except the case when fco 7^ and / 7^ 0. 



4.2. Reduction by two-dimensional subgroups. Despite the list of all twelve 
canonical forms of non-equivalent two-dimensional subalgebras of the ZK symmetry 
algebra given in Section I3.2.2( we shall only consider reductions of the equation 
by all pairs of distinct generators of the symmetry algebra (except those pairs 
of the form {^Fn^F-i} defined by the same type of generators). This is because 
space limitations precludes the treatment of all twelve cases in this paper. Also, as 
indicated in [11| , these reductions do not yield additional solutions when compared 
with those obtained from reductions by one-dimensional subgroups. However, they 
give rise to equations which are normally easier to solve. Each of the pairs of 
generators that we shall consider generates a subalgebra whose corresponding action 
turns out to have orbits of dimension two, and thus reduces the equation to an 
ODE. We shall always indicate to which of the classified canonical forms of two- 
dimensional subalgebras they correspond. 

The similarity variables in terms of which the reduced equations are expressed 
will be found by solving a system of first order partial differential equation of the 
form 

I . r — n 

(33) 

where each v and w are operators depending on the functions Fi and F2, respec- 
tively. Because of the arbitrary functions appearing in (|33p . it can be difficult to 
solve this system directly in certain cases. However, Equation (j33p becomes more 
tractable when it is transformed into the equivalent adjoint system of total differ- 
ential equation. The process for this transformation and the method for solving 
the resulting system of total differential equations are described in [501 HI] • 

One important aspect of Equation ([33)) is its integrability condition. If we con- 
sider for example this system with Vi?^ = Xg, and wj^j = (see Equation ([7])) 
and f{t) = 1, then we have Xg = (0,3,0, —5) and z/ = (1,0,0,0). We may thus 
ignore the variable y, in which case a function F[x, u) is an invariant if and only if 
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Xg • F{x,u) =^ 0. This in turn is equivalent to 

g' FJx.u) 

— = — r = a, wtiere a is a constant. 

9 K(x,u) 

This shows that the function g of the operator Xg most satisfy g ~ Ce°*, for some 
constant C. More generally, for an arbitrary pair {g, /} , the compatibility condition 
for Equation pS]) is given by 

"3/5'^-5V" + 3<?(/V + /5") = 0. 



This means that for any solution of to exist, we ought to have in this case 

K 



9 



/exp 



/ 



■dt 



(34) 



When the integrability condition is satisfied. Equation (j33p will yield two invariants 
w, and z, and in order to obtain the reduced equation we shall set w = w{z). 

Reduction by {vo,Xg} 

The corresponding Lie algebra belongs to type £2,9- The reduction formula is 



■u^9'v w{z) 2 

9y 



and the reduced equation is 



-2wg' 



This yields the group-invariant solution 



-V'9' 



9" 



.y'^9 '^99'. 

Reduction by {vo,y,J 

This Lie algebra belongs to the family of type £2,11 Lie algebras. The reduction 



formula is 



2x j/2 
17 ~ 2h 



w{z) 



2h 



The reduced equation is a generalized Riccati equation of the form 

w' + a(t)w2 + /3(i)w + 7(i) = (35) 

where a{t), P(t), and 7(t) are some functions of time. 
Reduction by {xg,y/i} 

This is a Lie algebra of type £2,4- The reduction formula is 



1 

2^ 



w{z)-2hg'x+{g'h' ~gh")!^ 



z^t. 



The reduced equation corresponds in this case to a mere integrability condition of 
the form 



~g'h' + 2hg" — gh'' 



0. 
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A first integral for tliis last equation is 



K = constant, 



9 9 

and this leads to a functional relation between h and g of the form 

r 

h = ag — Kg / g~^ dt, a = constant. 

The corresponding exact solution of ^ is thus 

^ w{t) ~ gH2xg' + yV)(a - K J g-^dt) 
" 2g^{a~ K J g-^dt) 

where w is an arbitrary function of time. For K = 0, it reduces to 

_ w{t)~ag^{2xg' + y^g") 
" " 2ag^ 

Reduction by {vo,z/} 

This is a particular case of type C2.1 Lie algebras. The equivalent system of 
differential equations takes in this case the form 
,,3 -f// 



&xyf-jir y , 

dy — di H dx 

^ 12xf 2x 



du — 



1 



12a;/ 

This leads to the reduction formula 



I2uxf' - Ax^f" - 2uy^f" - 2xy^f 



.,2 f (3) 



, u 
dt -\ dx. 

X 



u = 



y 



2x 



/' + /"/ 



6/2 

y2 6 ' 
and the reduced equation is 

12w - Uzw' + w'^ + {w + 2Az^)w" = 0. 
The solution for ui = is 



y 

6/2 



du 



-Augf + 2xf'g' - 2xgf" + y^g'f" - y^gf^^^ 



Reduction by {xg,z/} 

This Lie algebra belongs to type £2,10- The equivalent system is 

dy = 

1 

The integrability condition is 

-3/<?'2-5Y' + 35(/V + /5") = 0. 
Integrating this last equation once yields 

3/— - /' = K ^ constant. 
9 



dt dx. 

9 
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We must therefore have 



= /exp [ I ^dt 



f 



as we aheady indicated earher in 
The reduction formula is 



1 



w{z) 



x{K + /' 



y 



3J1/3 6/4/3 



J2/3 

and the reduced equation is the linear equation 



3 3 



with solution 



9z4/V + 6zi/V + — = 0, 
9 



(z) = -li^VAV3zi/V + 
18 



0,2 



for some constants of integration ai and 02. Substituting this expression for 
in the reduction formula yields 

_ 1 
18P 
1 

18/2 

Reduction by {y;i,z/} 

This Lie algebra is of type £2,12- The equivalent system is 



+ + 2/') - 3/(2i^x + 2xf + y V")) 



dy 



du = 



2yf hi2xf' + yyr 



2h , 

at H — ax 

yh' 



3/ 3y//i' 
h"i2xf + y^r) , (- W - 2x/" - 



6//1' 



dt — TT ax 
n' 



and the integrability condition is given by 

3//i'2 + 2h^f" - 3h{f'h' + fh") = 0. 
Integrating this equation once gives 

2f' -K _ h' 
3/ "77' 

where if is a constant of integration. Thus we must have 



h^C 



f exp I - / ^^dt 



1/3 



The reduction formula is 
1 



J2/Z 



w{z) 



y^f'{K-2n -6Kx + 6xf' + 3y^f 



2 fH-\ 



18/4/3 



18/1/3 



12xf + y'{K-2f) 

f4/3 



These yield the reduced equation 

+ 90Kw' + 1296w'2 + 36(36i« + Kz)w" = 0. 
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For K = 0, this equation reduces to 

w'^ + ww" = (46) 

with solution 



w — f3\/2z — a, {a = constant, f3 = constant). (47) 
The corresponding solution of ([T|) is 



18/2 ' "■'y—-' • ^ J ' • —^'-i V- J4/3 J- (4^) 

All the exact solutions to Equation ([1]) that we have found depend typically on zero 
arbitrary functions of time, and on at most two such functions and some arbitrary 
constants. Whenever the solution to a reduced equation was known, we were always 
able to readily obtain the solution to the original equation by a mere substitution 
into a reduction formula of the form (fT8|) . This situation is very different in the case 
of the direct method of Clarkson and Kruskal that we investigate in more details 
in the next section. Clearly, not all the group-invariant solutions are physically 
relevant for the sound wave propagation problem. However, they are certainly 
relevant for other problems modeled by the same differential equations when the 
boundary conditions change. 

5. Comparison with Clarkson and Kruskal's method 

By the direct reduction method of Clarkson and Kruskal [8] , all similarity solu- 
tions of the ZK equation of the form 

uit^x^y) = U{t,x,y,w{z)), z^z{t,x,y) (49) 

where f/ is a function of the indicated variables, and w{z) satisfies an ODE, may 
be found by substituting Equation ((49l) into Equation It is also argued in [7] 
that by a result of Clarkson and Kruskal, in order to find all such solutions, it is 
sufficient to look for u in the form 

u{t,x,y) = a + f3w(z) (50) 

where a and f3 are functions of t, x, and y. The latter equality is precisely our 
Equation (fT8|l that gives the general form of all possible similarity solutions that 
we have found thus far, in terms of the similarity variables. By making use of the 
reduced anzatz z^ ^ 0, Zhang et al. T' showed in this way that the most general 
ODE satisfied by w{z) has the form 



w'"^ + ww" + {Az + B)w' + 2Aw = -{Az + Bf (51) 



1 

3 

and raised the question of whether there is any connection between the direct 
method and the Lie classical method applied in the preceding section. 

One such connection could be determined by a way of mapping the reduced 
equation obtained by one method, to that obtained by the other method, taking into 
account that the direct method gives in principle the most general equation. Such 
a correspondence is not obvious because of the various changes of variables through 
which Equation ^ was obtained, except perhaps if equations are replaced by their 
equivalence classes under equivalence transformations. However, due to a well- 
known result asserting that two equivalent differential equations have isomorphic 
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symmetry groups [22j . these equivalences become easier to establish. The symmetry 
algebra of each of the reduced odes found in section l4?2l has dimension at least two, 
except for Equation (|38p whose symmetry algebra is generated by Vi = zdz + 2wdu, 
alone. For A 0, Equation (jSTj) is mapped after the change of independent variable 
y ~ z + B/A to the equation 

+ ww" + {Ay)w' + 2 Aw - \{Ayf = (52) 

which has the same symmetry algebra as Equation (|38p . The necessary condition 

for an equivalence between these two equations is thus satisfied. However, we've 

found that no linear fractional transformation of the form 

aiw + a2 

H= ■ , ^^a^x + ae, 

a^w + 04 

where ai, . . . , og are arbitrary constants with aia^^ — 0203 7^ 0, maps Equation (j52p 
to Equation ((38)) . This does not necessarily precludes the two equations from being 
equivalent under other types of transformations. 
For yl = 0, Equation ([5T|l reduces to 

- — + Bw' + w'^ + ww" = (53) 
o 

and the corresponding symmetry algebra has generators 'Vi — d^, V2 = zdz + wd^^. 
This symmetry algebra has the same dimension as that for the reduced Equa- 
tion ([l5|) from section 14.21 Equation (|^5)) with K — is exactly the same as 
Equation ([55]) with B — 0. For K ^ 0, the change of variables 



V — K'm(36w + kz), y — —Kz/36, with m — — 

36288 

reduces Equation (|45|l to an equation of the form 

B^ 

— —+Bv' + v'^ + vv" = 

with B = (—9 + 5V2T)-fsr/672. This shows the equivalence of the equations ([53]) 
and ([45]) . 

We thus see that by using the symmetry properties of the reduced equations, 
we can always map each of the symmetry-reduced equations to a sub-equation 
of ([5T|l . It should be noted that equations obtained by the direct method of Clarkson 
and Kruskal tend to be broad in nature and hence more difficult to solve. In 
addition, no symmetry or other properties of the solutions these equations are 
known, and there is also no guarantee that the solutions to the original equation 
is recoverable from those of the reduced equation, because the process of reverting 
back the solution generally involves solving other PDEs. When this process does 
give rise to an explicit solution of the original equation, it usually involves relatively 
complicated expressions. For example, the full expression of the exact solution 
generated by the reduced Equation (|46|) and given by ([47]) is easily found in Lie's 
approach and depends only on one arbitrary function of time and two arbitrary 
constants. The same solution is obtained by the direct method in [7| only through 
the solving a number of intermediary PDEs, and depends on no less than four 
arbitrary functions of time plus the two arbitrary constants. 

As for the reduction to (1 -I- l)-dimensional equations, the similarity solutions 
and corresponding reduced equations found in Section 14.11 are much simpler and 
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complete than those obtained using the direct method of Clarkson and Kruskal in 
[7], where both the equations and corresponding solutions depend on up to four 
arbitrary functions and are also usually defined only implicitly in terms of the 
solutions of some pdes. For instance, if a general reduced equation of the form |32)) 
were to be determined by the direct method, it would more likely be determined only 
implicitly and be expressed in terms of unnecessary arbitrary functions. Although 
the direct method also gives in principle the similarity solution as well as the reduced 
equation, its algorithm is relatively complicated and yields most often complicated 
and incomplete results. 

In the classical Lie reduction approach, the symmetry properties of each similar- 
ity solution is known and in some cases, especially in the case of solvable symmetry 
algebras, the solution to the original equation can be recovered by quadratures. 

It should however also be noted that even with Lie's classical method, having to 
revert back to the solution of the original equation by quadratures alone does not 
guarantee the explicit determination of this solution. Indeed, these quadratures 
may involve transcendental or integral equations that aren't easy to solve. This 
fact is common in the determination of group-invariant solutions and we give an- 
other example here by attempting to solve Equation (|45|) . for which the symmetry 
generators are given by 

Vl ~ {—'S6/k)dz + dw, and V2 = zdz + wdw 

These vector fields clearly generate a solvable Lie algebra with commutation rela- 
tions [vi, V2] = Vl. In terms of the rectifying coordinates 

y~2>Qw + Kz and r = — (Xx) /36 (54) 

for Vl, Equation (|^5)) reduces to 



-H + ISi?^ + 648i/^ + yH' = (55) 

where 

H=^. (56) 
dy 

Equation (|55p retains the symmetry V2, which is given by V2 = ydy + rdr in 
terms of the variables r and y, and by 

" (l-Hl8i7)i/3(36i/- 1)2/3 
in terms of H and y. In terms of the rectifying coordinates 

e - H (58) 

S^il±}^Ilf!lpILl^y (59) 

for V2, Equation ([55)1 reduces to 

e^(l + 180^/^(36^ - if^S' = 

with solution S'(^) — a, where a is a constant of integration. It thus follows from 
Equation (jSO]) that the solution H — F{y) of ((55)) is given implicitly by the poly- 
nomial equation 

aH^ -y^ + 5AHy^ - 23328H^y^ = 0. (60) 
Assuming that the function F{y) is known, it readily follows from Equation (|56p 
that 
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r = J F{y)dv - G(y) (61) 

for a certain function G(j/), and the substitution of (I54|) and (|44l) into this last 
equation leads to the solution to the original equation by an additional quadrature. 
The problem here with these quadratures is that the integral r = J F{y)dy is not 
easy to solve because of the complicated form of the function F. Indeed, one of the 
simplest roots H — F{y) of the polynomial Equation ((60|) is given by 

_ -S6ay^ + 839808/ + 2^^^{y^X^ + y/ay^X^)^/^ 
~ 22/3X(y3X2 + 7^7x3)1/3 

where X = a — 23328y^, and so the integral equation r = J F{y)dy does not appear 
to be obvious to solve, unless perhaps an appropriate change of variable to be found 
was performed. In reality, the difficulties with these quadratures are normally to 
be associated with the equation itself and not with Lie's reduction method. 

6. Concluding remarks 

In this paper we have given a classification of low-dimensional subalgebras of the 
ZK symmetry algebra into conjugacy classes and determined similarity reductions 
of the ZK equation to (1 + l)-dimensional equations and to odes. We have thus 
derived a number of new exact solutions to this equation in this way. Consequently, 
we have been able to compare the results obtained using Lie's classical method with 
those obtained in [7] with the direct method of Clarkson and Kruskal, and found in 
particular equivalence transformations between the reduced equations obtained by 
this two methods. This investigation shows that not only Lie's algorithm is simpler 
and much richer in properties, but it also yields simpler and more complete results 
as opposed to the direct method where both solutions and reduced equations are 
most often determined only implicitly by complicated expressions. 

The example of reduction by direct case analysis treated in Section 14.1.21 sug- 
gests that equations obtained by the direct method can also be viewed as reduced 
equations obtained by the Lie classical method by performing a reduction with a 
system of generic generators corresponding to Lie algebras having orbits of fixed 
dimensions, rather than performing the reduction with some specific generators. 
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